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Abstract 

Let G{x) be a C° function such that |G( 2 :)| < for |a;| < c, 

for constants K,c > 0. We consider spherically symmetric solutions of 
Og(j) = G{(j}) where g is a Schwarzschild or more generally a Reissner- 
Nordstrom metric, and such that (j) and are compactly supported on 
a complete Cauchy surface. It is proven that for p > 4, such solutions do 
not blow up in the domain of outer communications, provided the initial 
data are small. Moreover, |0| < G(max{u, 1})“^, where v denotes an 
Eddington-Finkelstein advanced time coordinate. 

The interaction between the geometry of black hole backgrounds and the be¬ 
haviour (in the large) of linear and non-linear waves plays a fundamental role 
in general relativity. The very stability of the simplest black hole solutions 
like Schwarzschild and Kerr could depend on subtle features of this interaction. 
Moreover, fine aspects of the behaviour of these waves, such as the phenomenon 
of decaying tails, have physical importance, both from the point of view of far 
away observers, and for those brave enough to cross the “event horizon”. 

In [^, we studied the problem of collapse of a spherically symmetric self- 
gravitating scalar field, i.e. solutions of the system 
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arising from spherically symmetric initial data, in the case where it is known a 
priori that a black hole formsd The main result of |S] was a proof of Price’s 
law, heuristically derived in EOlEl. The mathematical content of this “law” is 
a set of decay rates for the scalar field in appropriately defined null coordinates 
on the so-called domain of outer communications or exterior of the black hole. 
In particular, we showed decay |^| < in an Eddington-Finkelstein-like 

advanced time u-coordinate along the event horizon, and decay \r(j)\ < u~‘^ in 
retarded time u along null infinity. 

The problem of the self-gravitating scalar field studied in is non-linear. In 
particular, the black hole geometry depends on the solution, and even properly 
identifying the black hole exterior is a non-trivial part of the problem. The 
results of |H] can be specialized, however, to the easier problem of the study of 
spherically symmetric solutions to the equation 


□<,(/> = 0 


on a fixed Reissner-Nordstrom background spacetime These spacetimes 

constitute a two-parameter family of solutions to the Einstein-Maxwell equa¬ 
tions, with parameters mass M and charge e. The one-parameter subfamily of 
vacuum solutions defined by e = 0 is known as the Schwarzschild family. For 
parameter values 0 < |e| < M, these spacetimes contain black and white holes, 
and the so-called exterior region can be covered by coordinates (r, t) ranging in 
(M -I- 'fM'^ — e^, oo) X (—00,00) such that the metric takes the form: 


9 = - 




-1 

dr^ -I- ^2 ^ 


where da §,2 denotes the standard metric on the unit two-sphere. In decay 
rates were obtained both in the (static) t coordinate, as well as in Eddington- 
Finkelstein retarded and advanced time coordinates u and v. (The latter are 
indispensible for understanding the decay of the radiation tail along null infinity 
and the behaviour along the event horizon.) Previously, decay rates had been 
shown only in the static t coordinate mni under the additional restrictive 
assumption that the solution is supported away from the event horizon’s sphere 
of bifurcation.^ In contrast to the situation in Minkowski space, there is no 
decay rate in t which is uniform in r, in view of the geometry of the horizon. 

In this paper, we shall explore the implications of the decay proven in for 
the behaviour of small data solutions to the semi-linear wave equation 


□g<(. = G{<P) (4) 

on a fixed Reissner-Nordstrom background where the nonlinearity is 

assumed to decay sufficiently fast at 0: 

|G'(x)| < K\xf for |x| < c (5) 

^For the definition of the notion of black hole, see 1^. Black holes form for generic large 
initial data in view of Christodoulou’s proof of weak cosmic censorship [5| for spherically 
symmetric solutions to Q OJ. In the presence of an additional non-trivial Maxwell field, for 
which the results of also apply, the existence of a black hole follows from IHI . 

^We also note earlier results showing boundedness m and decay without a rate |22| . 
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for some constants K^c > 0. The study of small-data solutions to 0 under 
assumption © is often referred to as the “John problem”. The point of is 
that no special structure is assumed of the non-linearity; in particular, it is not 
assumed that G{x) = V'{x) for a nonnegative potential 1^ > 0. Our main result 
(Theorem^ states that for p > 4, small-data spherically symmetric solutions 
of 0) remain regular on the closure of the exterior region J~ (X+) n J+ {I~) 
of Schwarzschild or Reissner Nordstrom.^ Moreover, |(/)| < Cv~^ in advanced 
Eddington-Finkelstein time v, and, for r > ro > M -|- V— e^, |^|(r, t) < 
C{ro)t~^, but with C —*■ oo as ro —> M -|- ^/M'^ — e^! Note that our assumptions 
on the initial data in Theorem ^ are the geometrically appropriate ones; no 
special conditions are imposed at the sphere of bifurcation H\ n of the 
event horizon. 

The global study of the equation 0 on a Schwarzschild background was 
initiated in 12 EH!, in the case where the non-linearity G is given precisely by 
G{x) = x\x\'^. Although this is a special case of the phenomenology of this 
particular equation is completely different from the general case, as this non¬ 
linearity is of the form V'{x) for a potential R > 0. In view of the nonnegativity 
of R, the specific form of G,^ and the Sobolev embedding ^ L®, global 
regularity for large data solutions can be proven using only the energy estimate 
arising from contraction of the energy momentum tensor of (j) with the Killing 
vector dt-^ In particular, no appeal need be made to decay for the linearized 
problem. Under spherical symmetry, similar results can in fact be proven in the 
self-gravitating case, i.e. for the system 

- ^9tiuR = ( 6 ) 

(7) 

Rfiv 0,A (S) 

for a wide class of potentials. Specifically, if the potential is bounded below, 
i.e. V > —G > —oo, then a weak form of stability® has been proven in [ 7 ] 
for arbitrary spherically symmetric solutions of ( 0 - 0 , in the case were it is 
assumed a priori that the past of infinity has a nontrivial complement. If U > 0, 
the results of | 7 | imply in particular that null infinity is complete These 

®In stating results, we apply here standard notations from global Lorentzian geometry. See 
Hawking and Ellis in. The reader unfamiliar with these global concepts should note that 
our result implies in particular “global existence” with respect to Schwarzschild or Reissner- 
Nordstrom time t. 

"^in particular, the sufficiently low power 

^These methods have recently been extended 1191 to the case of a Kerr background, where 
dt is no longer everywhere timelike in the domain of outer communications. In this case, 
global regularity is proven even though the solution is allowed to grow. 

®i.e., that the solution does not blow up in the past of infinity 

^In the spherically symmetric self-gravitating case, an energy estimate is provided via the 
coupling by considering evolution equations for the Hawking mass m. 
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questions received attention [mini recently in the string theory community in 
the context of possible cosmic censorship violation. 

In contrast to the above, for the case of the “John problem” considered here, 
where we do not assume a special structure for G, results are only expected for 
small data solutions, and for p sufficiently large. We recall the situation in 
Minkowski space: The study of small amplitude solutions of the equation 

□</. = \^\P (9) 

on was initiated by F. John d, who showed that for p > 1 + v^, solutions 
to @ with sufficiently small, compactly supported initial data remain regular 
for all time, and moreover, decay to 0 along all inextendible causal geodesics. 
On the other hand, for 1 < p < 1 + F. John established the existence of 
small data spherically symmetric solutions which blow up in finite time.® This 
blow-up result was later extended by J. Schaeffer m to p = 1 J- \/2. 

Recently, it was proven in |3] that for equation on a Schwarzschild back¬ 
ground® for 1 < p < 1 J- there exist solutions arising from arbitarily small 
compactly supported data which blow up on the exterior region H 

In view of the restriction p > 4 in our Theorem^ it would be interest¬ 
ing to understand the behaviour in the interval 1 J- < p < 4, in particular, 
to determine the critical power of p separating stability and blow up.^® 

Besides addressing what is by now a classical problem in non-linear wave 
equations, our motivation for studying this “John problem” in the black hole 
context arises as follows: Experience from non-linear stability problems for the 
Einstein vacuum equations without symmetry strongly suggests that using de¬ 
cay rates of linearized fields will be an essential part of any argument. In the 
black hole context, implementing this approach introduces a number of tech¬ 
nical challenges, as the notion of “time” with respect to which decay is to be 
measured must necessarily be substantially altered. In this context, the present 
problem appears to concern the simplest geometrically well-defined non-linear 
equation where “global regularity” is indeed shown via decay for the linearized 
problem. Moreover, we believe that the estimates proven here clarify the role 
of the event horizon-and the celebrated red-shift effect-for ensuring stability. 
Indeed, it is not clear whether the results shown here could be obtained using 
only the (non-uniform) t-decay shown in mni, which do not precisely express 
behaviour up to the horizon. Further clarification of these issues may shed light 
on the problem of the non-linear stability of the Kerr vacuum solution. 

® Consideration of the general n-dimensional problem leads to the critical power pc defined 
as the positive root of the quadratic equation (n — l)Pc — (n. -|- l)pc — 2 = 0; this power first 
arose in work of W. Strauss. For analogues of F. John’s result for see irni and references 

therein. 

^The spherically symmetric geometric wave operator on Schwarzschild can be reinterpreted, 
after multiplication of (p by r, as the R^"*"^ wave operator, with an additional positive potential, 
and this is the approach taken in 1^. For extensions of the John problem to large classes of 
linear potentials, not necessarily positive, see |1 f)| . and references therein. 

^^For simplicity, we appeal in this paper only to the uniform u-decay estimate-not the 
complete set of decay rates-proven in [H|- We expect that a more detailed study should lower 
the constraint p > 4. It is not clear, however, whether the estimates of alone will be able 
to retrieve the result for p > 1 + y/2; for comparison with the case of see d 


4 



1 Schwarzschild and Reissner-Nordstrdm 

Let (A4,g) denote a non-extremal Reissner-Nordstrom solution corresponding 
to parameters M, e, with 0 < |e| < M. (We refer the reader to standard ref¬ 
erences, for instance m-) By convention here, the term Reissner-Nordstrom 
solution will always refer to a globally hyperbolic spacetime which is the Cauchy 
development of an asymptotically flat hypersurface S with two ends^^. Consid- 





r{q) = v' Area{q)/A'K, 


( 10 ) 


we have that 


g = g + r^dcr §2 


where dag 2 denotes the standard metric on the unit sphere, and where g defines 
a 1 -|- 1-dimensional Lorentzian metric on Q. The Lorentzian manifold {Q,g) 
can be covered by a bounded global null coordinate system (u, 1;). The range of 
such a coordinate system is depicted, as a subset of 2-dimensional Minkowski 
space, in the Schwarzschild (M > 0, e = 0) case: 


B 



11 


i.e., in the e ^ 0 case, we do not mean an analytic (or otherwise) extension thereof 
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and the charged Reissner-Nordstrom (M > |e| > 0) case: 


cnV '\C7i:+ 



These depictions are commonly known as Penrose diagrams. The curve S de¬ 
picts the projection to Q of a particular choice of complete Cauchy surface 
S CTW. 

We define the Hawking mass function m : Q —*■ K by 

^ (1-ff(Vr,Vr)), 


and the mass ratio 



The mass m is related to the parameter M by 


M = m + 

2r 

The sets J~{X^)r] and J”(Ig) n J+(Ig), the two so-called domains of 

outer communications, can each be covered by a coordinate system (r, t) with 
coordinate range (r_|_,oo) x (— 00 , 00 ), where r is defined by (llOll . and in which 
the metric g of Q takes the form 


g 




dV -I- 1 


2M 

r 



dr"^. 


Here r+ is defined by 

r+ = M +\/M^ - e2. (11) 

One immediately sees that the vector field dt is timelike Killing in these domains. 
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To define a null coordinate system, we first introduce the so-called Regge- 
Wheeler coordinate 

= r + r^(r+ — r_)“^ log |r — r+| — r?.(r+ — r_)“^ log \r — r_|, (12) 


where 

1 

CN 

II 

1 

This coordinate satisfies 

dr* 



We now set 

tt = t — r*, 


V = t -I- r*. 


With respect to these coordinates, the metric g takes the form 


g = — — ^)dudv = — ( 1-h 


dudv 


in J~ {l\)r\J^{I^) and J” ) n (X)) ); i.e. {u,v) are indeed null. Individu¬ 
ally, the coordinates u and v are each known as Eddington-Finkelstein retarded 
and advanced time coordinates, respectively.^^ Each of the domains of outer 
communications is covered by the range (— 00 , 00 ) x (— 00 , 00 ). 

The notion of asymptotic flatness is captured by the fact that 1 — /i —> 1 
as r —!■ 00 . To examine the geometry of the horizon, first, we note that the 
timelike Killing vector field dt extends to a null tangent Killing vector field on 
the event horizon Ti. = vanishing at the sphere of bifurcation 

Thus 1 — /r = 0 along 7i, and moreover m and r are constant. In terms of the 
parameters, we have 


m = 2r = 2r+ = 2(M -h \/M^ - e'^). 


This vanishing of 1 — /i relative to the Regge-Wheeler (t, r*) and null (u, v) 
coordinates is captured by: 


Cie 



{v — u) 


Cie 




< C'i(r-r+) 

< 1 - M < C' 2 (r - r+) < (726 ^ 


(13) 


(726 



{v — u) 


for small enough —00 < r*, and for some positive constants, Ci, Ci, C 2 , and 
C 2 , depending only on M and 6. 

^^The term Eddington-Finkelstein coordinate system, however, typically refers to the hy¬ 
brids (u,r) and (v,r). 
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2 Local existence 


Let G be a bounded C°° function of its argument, let {A4,g) be as in the 
previous section, and let E be a Cauchy surface. We have the following 

Proposition 1. Let (jjQ and (pi be C°° functions on E. There exists a unique 
globally hyperbolic subset V G M, and a C°° function p ■.'D such that 

1. G D is a Cauchy surface, and P\y: = po, N°‘pa = pi, where N is the 
future-directed unit normal to E. 

2. OgP = GiP) 

3. The set T> is maximal, i.e. ifUGL) is globally hyperbolic, with ^ : I? —> R 
satisfying the above, then D = 7). 

We shall call D the maximal domain of existence of the solution. Moreover, if 
E is spherically symmetric, and po and pi are constant on the SO{3) orbits, 
then D is spherically symmetric and p is constant on the SO{i) orbits, and thus 
descends to a function on Vq = V/SOlS). 

In the spherically symmetric case, the regularity assumption in the statement 
can be lowered to G € G;°^, po S Cl^^, pi S G;°^. The resulting solution p is 
then G^. In a slight abuse of notation, we will denote Vq again as V. It turns 
out that blowup can be easily characterized in terms of the L°° norm of p, 
i.e. we have 

Proposition 2. Let q G QDV. If J~ [q) GV %, and 

sup \P{x)\ < oo, 

x^j- ((j)n'D 


then q G T>. 

A similar statement can be made about the past. The proofs of the above 
proposition are completely standard and are thus omitted. 

3 The main theorem 

The main theorem of this paper is the following: 

Theorem 1. Let H be a spherically symmetric Cauchy surface in a Reissner- 
Nordstrom spacetime {Ai,g), with parameters 0 < e < M, and let po, pi be C^, 
G° functions, respectively, on E, constant on the 50(3) orbits, and supported in 
E n {r < R}, for some R < oo. Let T) denote the maximal domain of existence 
for the initial value problem defined by Proposition 0 where G is assumed to 
satisfy 0 with p > 4. Then there exists an e, depending only on M, e, E, and 
R, such that for 

I0o|ci < e, 

^^The norm above refers to the induced Riemannian metric on S. 



\(l>i\c° < e, 

T) contains the closures of the two domains of outer communications, 

(j-(i+) n J+{ij)) u n J+{Ib)) c V. 

Moreover, 

\(j){u,v)\ < C(max{t;, 1})“^, 

and, for each fixed tq > 

\4'{r,t)\<C{ro){\t\ + l)-^ r>ro, 

where v and t denote an Eddington-Finkelstein advanced time coordinate and 
a static time coordinate, respectively, on either of the two domains of outer 
communications. 


4 Proof of Theorem [TJ The bootstrap 

In view of simple arguments of Cauchy stability, and the fact that the data are 
compactly supported on S, it follows, that given any to, say t = 1, then for e 
sufficiently small, we have 

{t = l}cV 

and 

|</'|i=i|ci < C (14) 


\N(f\t=i\co < e. (15) 

Here again, the norms refer to the induced Riemannian metric of {t = 1}, and 
N refers to a unit future directed normal. The relation of 5 = T,/SO{3) and 
{t = 1} could be given, for instance, as follows: 



i 


Note 4)\t=i will not in general be compactly supported in (r+,oo). It will be 
supported in (r-+, R] for some R < oo. Moreover, we have e = e(e, S, R, M, e) —> 
0, as e ^ 0. 
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In proving our theorem, it clearly suffices to show that 


V D J+(s) n J-{iX) n J+{Ia). 

Define the set 6 C D n J“ {l\) n J+ ) by 

B = \q € V f] J~ {X\) n J^(X^) : sup v+{x)\<f>{x)\ < 5 > (16) 

[ x^j-{q)n{t>i}nv J 

where denotes max{u, 1}, and for a constant B to be determined later. 

This subset is clearly non-empty and open for B > €V+{R, 1). In view of 
Proposition [21 proving Theorem ^ reduces to showing that B is closed. This 
would follow immediately from 

Theorem 2. For q G B, we have 

sup v+{x)\(j){x)\ <—. 

xeJ-{q)n{t>i}n'D ^ 

To prove this, we shall view our non-linear problem in as a solution of 

= X, (17) 

where F is to satisfy the estimates arising from (P|) if F is set equal to the 
non-linearity. We thus turn in the next section to the study of this linear 
inhomogeneous equation. The proof of Theorem (21 will then follow easily in 
Section ini 

5 The linear problem 

First some formulas: Written with respect to Regge-Wheeler coordinates, equa¬ 
tion itTTl takes the form 

(1 - - r-^drAr^dr.-tp)) = F. (18) 

Alternatively, relative to the null coordinates [u = t — = t + r^) equation 

m can be written as a first order system 

dv{rduip) = ^K f, (19) 

Zr 4 

duirdy-ip) = ^——^rdutp + ^——^rF ( 20 ) 

2r 4 

for the pair of functions 0 = rdyt/j and by a single equation 

dxduiri;) = - (^1) 
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We also record the following form of equation crucial to understanding the 
so-called red-shift effect 


dy 


rduif 

1 - /i 


\ r ) I- n 


rdyip 

2r 


—rF. 

4 


( 22 ) 


This effect has been discussed in |S]. 

From m and TO-CT we can derive the following energy estimates. Given 
(ui, ui), for t < {vi ui)/2, define 


J t — U^ 


' t — Ux 


= 2 


f t — U\ 


{\rdui’{t, r*)p -I- \rdyii{t, r*)p) di\ 


For ti = (ui -f ui)/2 > ^2 > ts, we have the following energy identity: 


p2t2 — Ui p2t2 — Vi 

_l_ / \rdytp{ui,v)\'^ dv / \rdutp{u,vi)\'^ du 

J 2t‘i — u^ J 2t^ — ui 


1 ^t2 y-t+vi 

/ / {1 - ^)F{t,r^)dt'ip{t,r^)r'^ dr^dt. 

^ Jt^ J t-U-i 


Hs Jt—ui 

The estimate refers thus to the region depicted below: 


(ui,ui) 



We also have the following characteristic energy identity in the characteristic 
rectangle [u 2 ,ui] x [u 2 ,ui]: 


\rdu'f{u,vi)\'^du - 


\rdy'tp{ui,v)\'^ dv = 


\rduifiu,vo)\‘^du ■ 


\rdyip{u 2 , v)f dv 


+ 


(1 — ^)F[u, v){du4’{u, v) — dyijj{u, u)) dv du. 


The reader can refer to the diagram of Theorem 01 


5.1 The homogeneous problem 

The linear homogeneous problem, i.e. the case F = 0, was studied in |^, where 
pointwise decay estimates were proven. In this section, we will present these 


11 







estimates in a form suitable for understanding the inhomogeneous problem via 
Duhamel’s principle. 

First we have the following preliminary result: 


Theorem 3. Fix tq > r+, and let {ui,vi) be an arbitrary point. Let tg < 
(ui +vi)/2, and consider the lighter-shaded region 


£ = J {ri,ti)nJ+{{t = to})n{r>ro} 


depicted heloni^^: 



Let Ip be a solution to 

DgP; = 0 

on J~{ui,vi) n J'^{{t = to})- Let {u 2 , V 2 ) G £■ Then there exists a constant A, 
depending only on rg, e, and M, such that 


r\'tp{h,V2)\'^ < ^(?'|V’(u2,n2)P(23) 


r^\dy'ip{u,V2)\‘^ < A(^r^\dy'ipiu2,V2)\‘^ + , (24) 

for any point u > U2 with {u, V2) G £■ 

On the other hand, if{ug,Vg) is any point in J~ {ui,Vi)n J~^ {{t = tg}), then 


r\duip{u3,v)\ 

1- p. 


< A\ e '■^(“3.®) 


fy(5-«3) f r\duij{u3,V3)\ ^ [ 1 /:] 


for any v > V3 such that {u3,v) G J (ri,ti) n J'^{{t = to})- 
Finally, if 

C= sup {|r-!/'|(r,to), |r25«(r-!/')(r,to)|} 
£n{r>ro} 


(25) 


then 


rV’(r-,t)p + \r'^dyili(r,t)\^ + \r'^dy{rili)(r,t)\'^ < A , (26) 


throughout J (ci, ti) n {r > Tq} H {t > tg}. 

^"^Note that despite the diagram, we are not necessarily assuming that r(tti,ui) > rg. 
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Proof. Equations (E51) . lOl) . and iCTi can be easily deduced from the equations 
and the energy estimates of the previous section. For the estimate 
we write the homogeneous wave equation in the red-shift from given by 


rduif 

'l- 




1 — 2r 


Introducing the integrating factor e 


f:. 


^)d'’ 


and integrating in v we obtain 


rduHu3,i) ^ g-jj; ^{M-4)dv rdu^{u3,V3) ^ 

Since 

= (27) 

r r+ 2 

and r is increasing in the ?;-direction, we conclude 

r\dui’iu3,v)\ ^ r\du^{u3,V3)\ 

1 — ^ “ 1 — H 

f e ^r\dy'ilj{u 3 ,v)\dv. 

Jv3 

Applying Cauchy-Schwarz gives 


r\dutf{u3,v)\ 
1 - 


< g 2r^(. 


^iv-V3) r\dutlj{u3,V3)\ 


1 - M 


+ A[ |r3„■^/'(w3,^')P rfv) • 


□ 


We now turn to decay estimates. The full results of [S] are complicated to 
state, as different regions must be treated separately. In this paper we shall 
only use the following: 

Theorem 4. Let {ui,vi), (^ 2 ,^ 2 ) be points with {ui,vi) G J'^[u 2 ,V 2 ), i.e. with 
ui Pu 2 , vi >V 2 . Consider the characteristic rectangle 

TZ = J^{u 2 , V 2 ) n J~{ui,vi) 


depicted below: 



{U2,vi) 
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Let 'll) be a solution to □g'f/’ = 0 on 72.. Let 

C = max{ \rdy'il){u 2 ,v)\ + \r'^dy{r'il)){u 2 ,v)\, 

\{\- ^J)~^^u1|){u,V2)\] ■ 

Then there exists a constant A depending only on rg = r(u 2 ,V 2 ), e and M such 
that 

ll){ui,Vi) < AC{vi - V2)~^. 

The above theorem concerns a characteristic initial value problem. We can 
combine Theorems 13 and 01 to derive decay estimates for solutions of the homo¬ 
geneous Cauchy problem. 

Theorem 5. Let {ui,vi) be an arbitrary point, and and let tg < (ni -I- ui)/2. 
Let Ip be a solution to 

agP) = o 

on J~ (ui,vi) D {{t = tg}). Let 

C = sup{(l - yi)~^r‘^\dtil)\, (1 - p)^r‘^\dril)\,r'^\dv{r'ij))\,r\'il)\}. (28) 

t = to 

Pick a point (u 2 ,V 2 ) G J~{ui,vi) C {i = tg}: 



It follows that 


\'>Piui,vi)\ < AC{vi - V2) ^ 
where A depends only on rg = r{u2, V2), e and M. 

Proof. We first show that < AC^. Note that dr,ip = (1 — iTjdrip. 

2 


(29) 


E, 


(ui.vi) 

*0 


r — to+Vi 

I r^{\dtip\‘^+ \dr,ip\‘^)drt 

Iq — Ui 

r{r^— — to+vi) 

{l-yi)-^r^{\dtiP\^ + \dr,iP\^)dr 

r{r,=to — ui) 


< 


r-^ dr < C^rf^. 


fr+ 
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We apply now Theorem |3 From (I25II and 12811 . the estimate 


1(1-^) ^du’4){u,V2)\ < e '■^(“'" 2 ) '"”^“^^1(1-^) ^dui){u,vo{u))\+AC 


< g r-‘{-u,V2) 


(v2-Voiu)) 


{1- fi)-^C + AC 


(30) 


follows for U 2 < u < Ui, where Vg(u) denotes the z;-coordinate of the intersection 
of the outgoing constant-u ray with t = tg. Since, 


V 2 -Vg = 2(ro* - r*(u,U2)), 


we have by for sufficiently small r* = r*(u, V 2 ), r = r{u, U 2 ): 


< 


< 

< 


A{1- 

Ae 

(r-y2r-+)(r+y2r+)(r-+-r_) ^ 

A 2 ?’^ r 2 ^ 

Ae + 

j[g-r,(r-V2r+) 

A. 


(31) 


(Note that by convention here, A does not always represent the same constant.) 
Since the above estimate is also clearly true for large r, putting together pHl 
and m yields 

sup {1 - fj,)~^\du’tp{u,V 2 )\ < AC. (32) 

U2 

On the other hand, from OH), we have 

\dyi>{u 2 ,v)\ < ACr~^, 

\dy{rip){u 2 ,v)\ < ACr~'^, 

for vi < w < U 2 , in view of the fact that r > rg for such {ug, v). Our result now 
follows immediately from Theorem 0] □ 

Remark 1. Observe that the constant C in the above theorem is dominated by 
the (Riemannian) C^ norm ofr^ij), 

C < |rV|t=tolci- 

Remark 2. The quantity du'4^/dur is clearly independent of the choice of re¬ 
tarded time u and coincides with —(1 — fi)~^du4’ for Eddington-Finkelstein u. 
As one can choose a global regular null coordinate system, such that dyV < — 1 
along the ray u = U 2 in a neighborhood of TL\ O {u = M 2 }; one could have 
alternatively derived OH by appealing to Cauchy stability. 
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5.2 The inhomogeneous problem 

To prove Theorem |21 we must derive estimates for solutions of the equation 

= F, (33) 

with vanishing initial data at t = 1. Let ipir, t; s) denote the solution of = 0 
evaluated at (r, t) with initial condition 


'0(-,s)=O, dtipi-.s) = {l- n)F{-,s). (34) 

Using representation m we can then write the solution of as 

'I'(r, t) = J s)ds. (35) 

We prove the following theorem: 

Theorem 6. Let {ui,vi) be given, let dt satisfy on J~{ui,vi) H {t > 1}, 
with 'h = 0, 9(4' = 0 on J“(ui, ui) H {t = 1}, and let a > A. Then there exists 
a constant A, depending only on M, e, and a, such that 

|^'(Mi,ni)| < A(ui)“^ sup vf\F{u,v)\, (36) 


where v+ = max(n, 1). 

Proof. Let ri = r{ui,Vi), U = (ui + ni)/2. We will fix some vq > r+, to be 
determined later, and distinguish the cases ri > tq and ri < rp. 

5.2.1 Estimates in the region r > ro 

Consider first the case ri > cq. We have vi = ti + rl > + Tq. 

Let us take s > 2vf, for 0 < w < 1. To compute V'(ri,ti; s), we decompose 

J“(ri,ti) n {t = s} = {(uju) : u = 2s — V, vi — 2{ti — s) < v < ui} 

= /i U /2 


as depicted: 
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i.e. we have 


11 = {{u,v) : u = 2s — V, Vi — 2{ti — s) < v < Vi}, 

1 2 = {{ujv) : u = 2s — V, v'^ < V < vi\. 

Note that on the set Ii , 


2r* = V — u = 2{v — s) < —2vi 
and therefore, by we have 

(r^-r_) ^ 

1 - ^(r, s) < (726 \ (r, s)€/i. 

By the above estimate, (P|) and we have 

= A'^ f (1 -//)|F(r*, s)pr^ dr 
7/,u/2 


<71M e 


(r+-r_) ^ 

——’'i 


sup |F(r*, s)p + vf sup |F(r*, s)p , (37) 

h I 2 / 


where in the last inequality we used that r < C'vi on 12- 


On the other hand, for s < 2vf we can apply Theorem for the points 
(ri,ti) = (ui,ri) and (ro,s) = (u 2 ,U 2 ), where 


V 2 = s + rQ < 2vi. 

(In the last inequality, we assumed that rp is chosen in particular so that rp* < 
0.) Then, 

\i^{ri,ti;s)\ < Arf^sup(l - /r)5r^|F(r, s)| 

r 

< sup {1 — ^)^r^\F{r, s)\ + Av^^ sup (1 —/r)^r^|F(r, s)|. 

r*<-f '••>-§ 

For points r* < — estimate US gives 

_(r+-r_)^ 

1 - M < ^26 , 

while for r* > — we clearly have 

s 

r = r* + s > 
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Thus, for (5 > 0, 


fi; s)| < Av^^e ^ sup r^|F(r, s)| 

r.<-f 

+ Av^^s~^~^ sup (r, s)|. 

r.>-f 

Applying ll^ . and the estimates m and (EHli, we obtain 

- s 

|^'(?'i,ti)| <^iAe ^'’+ sup \F{r,s)\ + AtiVi sup 

J~iri,ti) J- (ri,ti)n{v>v'^} 

+ Av^^ sup s)| 

J~ (ri,*!) 

<(?;i)J sup |F(-«,u)| + (ui)+^ sup u++^r^|F(r, s)|, 


(38) 


|A(u,?;)| 


J (ri,ti) J (ri,ti) 

where we have used the fact that for ri > rg > r+, we have (fi)+ >~ ti, with 
constant depending on rg. 

Clearly, there exists a constant H, depending only on M and e, such that 
for r > H, we have v = t + > r in the region f > 1. It follows immediately 

that for all r, we have r < H + \ v\ in the region t > 1. We obtain, for any a > 0, 

\^{ri,ti)\ <A{vi)l~‘^°‘ v^\F{u,v)\ 

J~ (ri,ti)n{v>tJ^} 

+ A(z;i)^^ sup v^^\F{u,v)\ (39) 

<A{vi)^ ““ sup u“|F(m,u)| 

+ A{vi)^^ sup v^^\F{u,v)\. 


Choosing 


we obtain 


a = A + S, 


2(4 + 5) 

sup u“|F(m,u)|. 
J~ 


(40) 


5.2.2 Estimates in the region r < rg. 

Consider now the case r_|_ < ri < rg. In this region we shall use the red-shift 
effect manifest^® in the equation lUl): 


5, 


rd„At 1 


1 - M 


+ 1^-7 


rd,,'^ 


1 - M 


2r 


1 

—rF. 

4 


^®See 1^. The point, as we shall see, is that we can extract decay in v for du'^, even though 
we are integrating in v towards the future. 
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For fixed u, let Vt{u) be defined by t{u,Vt{u)) = 1. Using the integrating 
factor '> and the vanishing of the initial data on {t = 1}, we 

obtain 


rdu'i>{u,v) 

1- ^J. 


1 

4 



r =(M-eV ^)(2d,^i>{u,v')+rF{u,v')'^ 


dv' 


Integrating by parts we derive 


rdu'^iu,v) , 

+ f '■ '> ( — 2r~'^{M — e‘^r~^)'^{u,v') + rF{u,v') \ dv'. 

Jv,{u) ^ ' 

Since in the region under consideration, we have the bound < r < vq, and 
also, 


M - > M -,_ 

r M + VM2 - e2 


= - e2 > 0, 


it follows that 
|9„4'(u,i;)| 


A|4'(w,u)|+A [ e-<''-'''H\'i>{u,v')\ + \F{u,v')\)dv'. 

Jv,{u) ^ ' 


Given u, let u*(u) be defined by r{u^{v),v) = rg. Then 


|'I'(u,'y)| < |'I'(M*(r;),r;)| +A 


(1 — n{u', v))|'I'(u', u)| du' 


> (tj) 


+ A 





dv'{l — v)) du'. 


We now observe that in the region r_|_ < r < tq, we can estimate |1 — ^| < S for 
some 6 that can be made arbitrarily small provided that rg is sufficiently close 
to r_|_. Moreover, 


|1 — v)| du' = r{u, v) — rg < 5. 

' (tj) 

Therefore, for any u > Ut{v), 

|'I'('u,w)| < |'I'(it*(?;), r;)| + A5 sup 

u' G )'li] 


+ A5 sup 


u'^[U:t<{v),u] J V* {u') 




(41) 

dv' 
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Split 


Jv,{u') ^ ^ 


=-<v-v')( I / / / / 


u',v')\ + \F{u',v')\)dv' 


l-v/z 
^ v* (u' ^ 

+ £ ^;')| + \F{u', i;')|) dv' 


and estimate 

u/2 




/d/2 


< 


Ae sup ^|\E'(m','!;')| + |F(u',i;')|y 

v'^[v*{u'),v\ ^ ' 


+ \F{u\v')\^ dv' 

<A sup ('|d;'(M',u')| + |F(u',u')|)- 

f'' ^[v/2,v^ ^ ' 


V' ^\v/2,v 


Returning to (I41II . we derive 

|d/(it, u)| <|'I'(u,(u), u)| + A(5 sup d/(it',u) 

u'^ [li* ('l/),1i] 

+ A5 sup sup u')| + |F(u', u')|^ 

u'^[u:^{v),u] v" G[v/2,v] ^ ^ 

+ A5e“°’' sup sup ^|'I/(u', u')| + |F(m', u')|y 
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Thus, 


u)| u)| + sup 

u'^ [u^ (i;) ,u] 

+ v^AS sup sup |'I'(u',u')| 

■u'€(i)) v'^[v/2,v] 

+ Av^6e~™ sup sup |'h(u',u')| 

u' ^[u*{v) ,u] v' ^[V:t,{u'),v\ 

+ v^+A5 sup sup \F{u',v')\ 

u'G[u^{v),u] v' ^[v/2,v] 

+ Av^Se~™ sup sup |F(u',u')| 

u' ^[u*{v) ,u] v' ^[v^{u') ,v] 

u)| + sup 

u' (v) ,u] 

+ A5 sup sup z;')| 

u'G[u^{v),u] v'^[v/2,v\ 

+ AS sup sup |'I'(u', z;')! 

u'^[u^{v),u] v' ^[v*{u') ,v\ 

+ AS sup sup {v'_^_)^\F{u',v')\ 

u'^[u (v) ,u] v' ^[v/2,v\ 

+ AS sup sup |T'(m',z;')| 

li'G [li* ,li] v' ^[v^,{u') ,v\ 

<v^\'^{ut{v),v)\ + ASv^ sup 

u' G [tt* (v) ,u] 

+ 2AS sup sup z;')| 

F2AS sup sup (z;^)^|F(u', u')|. 

u' ^[u:i,{y),u] v' ^[v:i,{u'),v] 

The constant A is independent of the choice of tq and thus of S. A fortiori^ the 
above estimate applies also when \v^'i'{u, ?;)| on the left hand side is replaced by 
sup„' 6 [u 4 „),„] sup„,g[„^|( z;+)^4'(m',z;')I- Thus, if tq is chosen close enough 
to r+, so that S is sufficiently small, we have at our point (ri,ti) the estimate 

< 2(vi)^^f sup i;^|^'('u*(z;),z;)| 

'' vG [tJ* ,1^1] 

+ sup sup v!^IF(u, v)l), (42) 

uG[u^ (vi),uij vG[y* (u),vij ' 

for any /3 > 0. Here v* denotes the u-coordinate of the intersection of the curve 
r = rg and {t = 1}. Observe that the first term on the right hand side of 
the above inequality can be rewritten as sup]^<«j^ u^|'I'(ro, t)|; the supremum 
is thus taken over dz evaluated at points in a subset of the region r > rg. We 
choose 

P = l. 


21 


From 14011 and 14211 . we infer now that 

max i;“|F(m,v)|, (43) 

for all (ri, ti). □ 

6 Proof of Theorem [2] 

Recall the set B and the constant B defined by itTHll . Let (ri, ti) = (mi, iii) S B. 
Set F = G{(j)) and consider 'I' a solution of (IHHl in J~{ui,Vi) with vanishing 
data at t = 1, and let ip he a solution of the homogeneous equation with data 

"^1 J” 1} 4^\j~{ui^vi)n{t—l}: 

dt4’\j-(ui,vi)n{t=i} = dt'P\j-{ui,vi)n{t=i}- 

We have that 

(j) = tjj + (44) 

in J~{ui,vi) n J^{{t = 1}). 

By Theorems 121 and 01 it follows that 

|^(miWi)I < (45) 

where E = E{e, M, e, E, i?) ^ 0 as e —> 0. If 0 < i? < c, then by assumption ©I 

|G(0(u,u))| < K\(l){u,v)\^,\l{u,v) G B 

Thus, by the definition m of the set B, we have 

\F{u,v)\ < KBPvpP. 

Therefore for p > 4, choosing a = p and applying Theorem |21 we have 

|4'(wwi)l (46) 

and thus, by lEI) and (1^ . 

muuvi)\ < A{vi)p^KBP + E{vi)p\ (47) 

Let us define B = 4if, and let us require E (and thus e) to be sufficiently small 
so that AKBP < -j. Inequality (IT7I) then yields 

as required. 
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